Baryon Rapidity Loss in Relativistic Au+Au Collisions 
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An excitation function of proton rapidity distributions for 
different centralities is reported from AGS Experiment E917 
for Au+Au collisions at 6, 8, and 10.8 GeV/nucleon. The 
rapidity distributions from peripheral collisions have a valley 
at mid-rapidity which smoothly change to distributions that 
display a broad peak at mid-rapidity for central collisions. 
The mean rapidity loss increases with increasing beam en- 
ergy, whereas the fraction of protons consistent with isotropic 
emission from a stationary source at midrapidity decreases 
with increasing beam energy. The data suggest that the stop- 
ping is substantially less than complete at these energies. 



PACS number(s): 25.75.-q, 13.85.Ni, 21.65.+f 



Nuclear matter is believed to be compressed to high 
baryon density (ps) during central collisions of heavy 
nuclei at relativistic energies 0-0]. In the interaction 
region of the colliding nuclei, nucleons undergo collisions 
which reduce their original longitudinal momentum. This 
loss of rapidity is an important characteristic of the re- 
action mechanism, and is often referred to as stopping 
H . The rapidity loss for beam nucleons has been exten- 
sively studied in p+A reactions || and more recently in 
heavy- ion reactions [|6|— |To[| . Relativistic heavy-ion colli- 
sions are unique in the sense that secondary collisions of 
excited baryons are expected to contribute to the rapid- 
ity loss leading to simultaneous stopping of many nucle- 
ons within the interaction volume. The ps thus reached 
may be large enough to induce phase transitions, such 
as quark deconfinement and/or chiral symmetry restora- 
tion The observation of large numbers of baryons 
at mid-rapidity is indicative of compression to high pb, 
but the quantitative connection to ps is only possible 
via model calculations. An important key to our under- 
standing of this phenomenon is systematic data for the 



rapidity distributions of baryons over a broad range of 
conditions, such as beam energy and centrality. 

In this Letter, we present an excitation function of 
the centrality dependence of proton rapidity distributions 
from Au+Au collisions at 6, 8, and 10.8 GeV/nucleon. 
At these energies, the measured ratio of antiprotons to 
protons is very small (< 0.03%) Q and the production 
of protons from A decay contributes less than 5% to the 
total yield ]T^j. The measured protons can therefore be 
considered to directly reflect the distribution in rapidity 
of the initial baryons (assuming that the neutron rapidity 
distribution has a similar shape). The measured rapid- 
ity distributions at all three beam energies show a sim- 
ilar evolution with centrality. They are clearly bimodal 
in shape for peripheral collisions and change to shapes 
which for central collisions may still be bimodal in na- 
ture. This suggests that the degree of stopping in central 
collisions is not complete, even in this heavy system. In 
addition, it is shown that the protons do not end up 
as an isotropically emitting source at mid-rapidity, but 
retain a fair degree of their initial longitudinal motion. 
This phenomenon is represented by a parameterized fit 
to the data with a superposition of longitudinally moving 
sources. 

Experiment E917 measured Au+Au reactions at beam 
kinetic energies of 6, 8, and 10.8 GeV/nucleon at the 
Brookhaven AGS, and was the final experiment in the 
series E802/E859/E866/E917 @g|). The experimen- 
tal apparatus in E917 consisted of a series of beam-line 
detector arrays which were used for global event charac- 
terization, and a large rotatable magnetic spectrometer 
used to track and identify particles. The tracking sys- 
tem of the spectrometer consisted of a series of drift and 
multi-wire ionization chambers, which bracketed cither 
side of a dipole magnet, followed by a segmented time- 
of- flight wall of vertical scintillator slats. The data pre- 
sented here were taken with a trigger that required at 
least one track in the spectrometer. The centrality of an 
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event can be selected through either of two approximately 
equivalent methods: (1) from the multiplicity measured 
by a large acceptance device called the New Multiplicity 
Array [p.6| , or (2) from the energy deposited in the zero 
degree calorimeter jl6| . The data were sorted off-line into 
different centrality classes, where normalization was pro- 
vided by prescaled interaction triggers. For each beam 
energy, the five event classes are reported as the percent- 
age of the total interaction cross section (<7i n t = 6.8 b), 
corresponding to (0-5)%, (5-12)%, (12-23)%, (23-39)%, 
and (39-81)% ((39-77)% for 10.8 GeV/nucleon) of cr int . 

The systematic uncertainty on the normalization of 
the measured invariant spectra and rapidity distributions 
is dominated by the uncertainty of the single-track effi- 
ciency and the loss of tracks due to hit-blocking. The 
tracking uncertainty increases from 5% to 10% towards 
mid-rapidity. For peripheral collisions, there is a 10% 
uncertainty in the cross section, which decreases to 5% 
for central collisions. These uncertainties lead to a total 
systematic uncertainty of 15% independent of centrality, 
with a 5% relative uncertainty across beam energies. 
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FIG. 1. Invariant yields of protons as a function of trans- 
verse mass for ten rapidity intervals for each centrality class 
of Au+Au collisions at 8 GeV/nucleon. The most backward 
rapidity in each panel is plotted on the correct scale, while 
successive spectra have been divided by ten for clarity. The 
errors are statistical only. The curves are Boltzmann fits de- 
scribed in the text. 

The measured invariant yields of protons from Au+Au 
collisions at 8 GeV/nucleon are shown in ten rapidity 
intervals in Fig. | l| as a function of the transverse mass, 
m t = y/pt + tUq. The spectra at 8 GeV/nucleon are 
representative of the quality of the data at other beam 
energies. The errors are statistical only. The data were 
fit with a Boltzmann- form function in wit- 
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where T and dN/dy are free parameters. The data points 
were weighted according to their statistical errors. As 
has been previously reported, for Au+Au collisions at 
10.8 GeV/nucleon j|, the proton spectra cannot be sat- 
isfactorily described by a single exponential function. In 
contrast, the above form reproduces the spectra well with 
X 2 per d.o.f in the range 0.5 - 2.0, and provides the inverse 
slope parameter (T) and the rapidity density (dN/ dy) for 
each rapidity interval. 

The centrality dependence of the rapidity density at 
each beam energy is shown in Fig. |2[ The data are shown 
relative to the rapidity of the center-of-mass of the sys- 
tem (j/cm) which equals 1.35, 1.47, and 1.61 for beam 
kinetic energies of 6, 8, and 10.8 GeV/nucleon, respec- 
tively. The data points reflected about mid-rapidity are 
shown as open symbols. The errors were calculated from 
the fitting procedure. The present values of dN/dy for 
the most central event class at 10.8 GeV/nucleon are in 
good agreement with the previously reported values from 
the E866 Collaboration 
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FIG. 2. Proton rapidity distributions for all centrality 
classes at all three beam energies from Boltzmann fits to the 
invariant cross sections. The open symbols are the data re- 
flected about mid-rapidity. The errors are statistical only. 
The curves represent double Gaussian fits to the data (see 
text), the centroids of which are indicated by the arrows. 

For each beam energy, a common trend is observed in 
the evolution of the shape of the rapidity distributions 
as a function of centrality. For the most peripheral event 
class the distribution has a minimum value of dN/ dy w 6 
at mid-rapidity. This concave shape persists to the next 
most central event class, corresponding to a centrality cut 
of (23-39)%, consistent with the expectation that most of 
the participant protons reside at beam rapidities follow- 
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ing these relatively peripheral collisions. For more central 
collisions, the rapidity distributions become progressively 
flatter, and begin to develop a broad maximum at mid- 
rapidity for all beam energies. The trend of the data 
suggests, however, that for the most central collisions, 
the distribution does not evolve to a single peak cen- 
tered at mid-rapidity, but rather is consistent with two 
components each displaced from mid-rapidity, or with a 
set of sources spread throughout the rapidity range. To 
qualitatively emphasize this point, we have fit each of the 
measured distributions with two Gaussian peaks centered 
symmetrically about mid-rapidity. These are shown su- 
perimposed on the data for the three most central event 
classes at 8 GeV/nucleon in Fig. || and clearly display the 
trend. This result, somewhat surprisingly, implies that 
complete stopping is not achieved at AGS energies and 
that the longitudinal rapidity distribution is a result of 
transparency. 
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FIG. 3. Fits to the proton dN/dy distributions for the three 
most central event classes at 8 GeV/nucleon illustrating the 
evolution of shape with centrality. The solid curves are the 
sums of two individual Gaussians (dashed curves). 

To quantify these results, we calculate the mean ra- 
pidity loss from the fits for the most central collisions 
(Sy) = Ubearn ~ (u) where (y) refers to the Gaussians cen- 
tered at positive rapidity. These values are listed in Table 
I. To compare with previous results for the most central 
data at 10.8 GeV/nucleon, we have also calculated the 
mean rapidity loss for protons in the restricted rapidity 
range —1.11 < y — y cm < and find a value of 1.07±0.05, 
which is in good agreement with the value of 1.02±0.01 
obtained for the (0-4)% most central Au+Au collisions at 
the same beam energy reported by Videbaek and Hansen 

The mean rapidity losses determined from our data 
show a systematic increase with incident energy. The 
mean rapidity loss does not, however, fully characterize 
all the features of the final rapidity distribution. De- 
spite the fact that a considerable amount of the origi- 
nal longitudinal momentum of the incident particles has 



been lost, the situation is clearly far from one where they 
are completely stopped. Fig. [| shows the dN/dy distri- 
butions and inverse slopes for the (0-5)% most central 
event classes at all three energies plotted as solid points. 
The dashed curves in Fig. ^ show the expected distribu- 
tion for complete stopping - i.e. isotropic emission from 
a source at rest in the center-of-mass system emitting 
protons with a Boltzmann energy distribution, the effec- 
tive temperature of which is adjusted to reproduce the 
inverse slope of the transverse mass spectrum at mid- 
rapidity. The effective temperature thus accounts for the 
combined effects of radial expansion and temperature of 
this source. 




FIG. 4. Measured (solid points) and reflected (open 
circles) proton rapidity distributions and inverse slopes are 
shown as a function of rapidity y — y cm for beam energies of 
6, 8, and 10.8 GeV/nucleon. The arrows indicate the tar- 
get and beam rapidities. The dashed curves represent the 
expected distribution for isotropic emission from a thermal 
source at rest in the center-of-mass system (y — j/ cm ), whereas 
the solid curves correspond to an optimum fit to the data for 
a uniform distribution of sources within a range of rapidities 
(j/cm —Vb)<y< (ycm + yb) with a Gaussian profile of T e s{y) 
centered at jy cm (see text). The parameters obtained from 
a least squares fit to the data are yb = 0.990, 1.086, 1.166, 
T o ff = 0.253, 0.267, 0.279 GeV/c 2 , and cr T = 0.697, 0.762, 
0.809 for -Bbeam = 6, 8, 10.8 GeV/nucleon, respectively. 

Whereas this gives a good description of the rapidity 
dependence of the inverse slope parameters, it completely 
fails for the dN/dy distributions. In fact, only a fraction 
of the observed yields can be accounted for in such a sce- 
nario. The yields can, however, be adequately accounted 
for by emission from a continuum of isotropic sources, 
uniformly distributed over a rapidity range y cm ± yb- In 
order to account for the rapidity dependence of the in- 
verse slope parameter it is, however, necessary to intro- 
duce a Gaussian rapidity dependency of the effective tem- 
perature T cS (y) = Til exp(-(y - y cm ) 2 /2cr|). The solid 
curves shown in Fig. ^ were obtained by a four parameter 
simultaneous fit of?/;,; T° ff ; ot \ and Nq (dN/dy normaliza- 
tion) to the experimental dN/dy and inverse slope T(y) 
values. Based on the fit cW/dy-distributions, the total 
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number of protons within the range — 2 < y — y cm < 2 
(1-2 in;dy) was found to be 155, 164, and 159 for E hcmu 
— 6, 8, and 10.8 GeV/nucleon, respectively. Since the 
total initial number of protons (158) is approximately 
accounted for, this suggests that the extrapolation into 
the unmeasured region is quite reasonable. 

The fraction of the observed protons corresponding 
to complete stopping (fi SO ) was calculated as the ratio 
of the areas under the dashed and solid curves, respec- 
tively. The values of fi so at each energy are listed in 
Table I. The values of the absolute rapidity loss from the 
distributed source fits are also listed in Table I. They 
are in close agreement with the values from the Gaus- 
sian fits and show the same systematic increase with in- 
creasing incident energy. This result is different from 
observations in p+A collisions in which the rapidity loss 
was found to be independent of beam energy ||. This 
would suggest an increased role of secondary interac- 
tions which increase stopping in heavy-ion collisions. In 
contrast, the fractional rapidity loss, (Sy) /Sy max , where 
Sy m ax = Vbeam - Vera, is found to be essentially constant 
with beam energy from 6 to 158 GeV/nucleon, similar to 
the trend noted in Ref. M. 

In summary and conclusion, we have measured pro- 
ton rapidity distributions for Au+Au collisions at three 
energies as a function of collision centrality. These distri- 
butions show a consistent evolution with increasing en- 
ergy and centrality leading to maximal rapidity loss for 
the most central collisions and the highest energy. This 
result suggests the importance of secondary reactions in 
the heavy-ion case, contrary to the situation in p+A col- 
lisions. Nevertheless, the degree of stopping of the inci- 
dent baryons is far from complete and only a fraction of 
the observed protons can be accounted for by the emis- 
sion from a stopped isotropically emitting source. The 
remainder still possess considerable longitudinal momen- 
tum. It is not possible to say unequivocally from the 
present data alone whether or not this corresponds to a 
situation where the colliding nuclei are, to some extent, 
"transparent" , or a fully stopped and compressed system 
has re-expanded. However, the systematic behavior of 
the centrality dependence of the shapes of the rapidity 
distributions suggests, surprisingly, that the former is the 
case. 
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TABLE I. Comparison between different measures of stop- 
ping in central heavy-ion collisions. The absolute {Sy) and rel- 
ative (5y)/5y max rapidity losses obtained from double Gaus- 
sian fits to the diV/dy-distributions are given in column 2 and 
3, respectively. The absolute and relative rapidity losses ob- 
tained from the fitted solid curves in Fig. ^(a-c) are given in 
column 4 and 5, respectively. Upper limits on the isotropic 
fraction /j so of the total dN/dy distribution (solid curves in 
Fig. [|(a-c)) are given in column 6. The last row lists parame- 
ters obtained from an analysis of 158 AGeV Pb+Pb collisions 
at 0-5% centrality [§). 
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"Fraction of isotropically emitted proton (dashed curves in Fig. ^(a- 
c)) of total number of protons (solid curves of Fig. Q(a-c)). 



^Statistical errors only. 

c Obtained from a double Gaussian fit, dN/dy = 85[exp(— {y— y C m + 
1.19) 2 /1.6) + exp(— (y — y cm — 1.19) 2 /1.6)] to the experimental data 
|. 

d Obtained using T; so = 0.26 GeV/c 2 , which reproduces (pt) at 

!/-|/cm=0|. 
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